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Introduction 





Motivation 


Typical analytical treatments of frequency selective 
surfaces (FSS): 


— Full-wave analysis 


Floquet theory, spectral domain techniques 
Commercial software 

Accurate, but full numerical implementation 
Expensive for iterative optimization 

Affords limited intuition for design 


— Lumped element / transmission line / circuit models 


e Typically treat FSS as bulk structure 


All individual unit cell interactions are aggregated 


Simple models often only applicable to a narrow 
range of geometries 


FSS circuit models from literature: 




















LC resonator [1] 


Individual modal 


Multilayer FSS [2] admittances [3] 


LC lel 
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Introduction 





Motivation 


Goal: develop an element-level circuit model 
— Structure FSS circuit models from literature: 


e Each FSS element represented by individual 
lumped components 




















e Coupling between any two elements represented ` 
by separate lumped elements 
c 
— Canonical (applicable to wide variety of designs) 
B LC resonator [1] 

— Facilitate the design of more complex FSS structures 

© Multilayer Multilayer FSS [2]  8admittances [5] 

e Reconfigurable ur 1 H ` 

. . . -LCI [CI 

— Inverse problem (FSS optimization) " T n T «w 


eL iLa 


e Optimization typically performed with genetic | MR | 
algorithm (GA), particle swarm optimization (PSO), 
etc. using full-wave techniques [4, 5] 








e Simpler methodologies would be desirable 
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Introduction 





Conceptual Model 
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Circular aperture FSS 

self impedance elements (Z.) associated with each unit cell 

Add coupling impedance elements (T network, Z, and Z;) between unit cells 
Remove FSS geometry, leaving a lumped element circuit problem 


(a 
(b 
m (c 
(d 
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Introduction 


Outline 


= Lattice representation of FSS 

= N-Port network representation 

"= Lattice network analysis 

= Computation of lattice admittance elements 

= Extension to multi-element designs 

= Extension to multilayer designs 

Example implementation: rectangular aperture FSS 

x FSS synthesis through lattice model 

= Polynomial representation of lattice admittances 

= Formulation of optimization problem through lattice model 

= Optimization examples 
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N-Port Network Representation 

Lattice Network Analysis 

Computation of Lattice Admittance Elements 
Extension to Multi-Element Designs 
Extension to Multilayer Designs 


Examples: Rectangular Aperture FSS 


FSS LATTICE MODEL 
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FSS Lattice Model 





N-Port Network Representation 








D | Aperture-type FSS 
Ao 
l + 
MAN 
` ni Ma Arbitrarily shaped apertures in a 
DE rs Ay metallic ground plane 
AL Each aperture assigned a network 
Lx Ay SC? port with a scalar voltage V, 
tr E S Mah ty : 
- sa Y,: self admittance of each aperture 
| I Ay | W nn: mutual admittance between 
1 + y 
y B X aperture 1 and aperture n 
(a) à w 
D a | 7-network between ports accounts 
ay Ay — P for coupling between any two 
, m2 [1 
" D et apertures 
| x AP. Port voltages (I.e. aperture fields) 
Il + Ay . . . 
,6 9 excited bv input currents (i.e. 
(b) : incident fields) 
I 
; 
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FSS Lattice Model 


Lattice Network Analysis 
i Equivalent input impedance 
, cy A "Active" input impedance [6]: 


Vir 


TA 

A » .) 
N V. ` Se >Z ij Z ij 
P» jA I (r,) 


| Finite difference equation: 











| ” HY, H! YV (v,)- $ Y, [V(r.) | - LV (v) 2 1(x.) 
HS | 


— 
d- 
"e 
ey d) d 
d 
Fee 
— 
N + 
Ex] 


AY oolution in terms of lattice 
til e? Ay Green’s function [7, 8]: 





^ Cie - XX G(ni.r')/() 
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FSS Lattice Model 


Lattice Network Analysis 


. Lattice Green's function 





Ao 
` , KR 9 Ay Green's function equal to Z-matrix 
: DK Ai element: 
I, | ; 
a.m " Z, - G(r.r;) 
p E | oe . l 
BU & Ein, 1 y Eigenfunction expansion of Green's 
Al 4 function [9]: 
E Ay 
al 18 9 iert, 
$ ei G(r.r,)- 5. P 
s LS e m—0 m 
E 4 CY» I ` 
l u Eigenvalue problem: 
Hm 
(b) l LV, (v) - A,V, (r) 
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FSS Lattice Model 


Lattice Network Analysis 





Uniform Plane Wave Incidence 





Ged 
8 i Y . ° . 
^ D : = Under uniform incidence: 
| 7 | Ay — Equal voltages excited on all apertures 
"mE GZ: — Uniform input currents 
S = Lat 
Ka Ay À m = By symmetry: 
l + a Ay EN le + y 
y, Es UU vm — No current flows between ports 
= | Ay CF through pi-networks 
OE 
l bk Cr Ay — Lattice collapses to equivalent 1-port 
(a) : network with infinite admittance sum 
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FSS Lattice Model 


Computation of Lattice Admittance Elements 





Handling infinite summation of admittances 


Rewrite infinite sum of x o o 
admittances as two- |. Y^ -Y, «4» Y,-» » Y, 





dimensional series: ba a 
Corresponding partial sr wr 
YMN — » » P 
sums: = 
Central aperture (self admittance) 
Bl Level 1 coupling 
MI Level 2 coupling 
B Level 3 coupling olowly converging two-dimensional 


series (similar problem arises in antenna 
array theory [10]) 





XT | TEXAS A&M e 


I V ERSITY 





FSS Lattice Model 


Computation of Lattice Admittance Elements 





Handling infinite summation of admittances 


Real part of normalized TI” sums 


Solution: 2-dimensional Shanks transformation 


Transformation computed as a determinant 
comprising the partial sums of the series [11, 12]: 














A A, ta l éi l 
y'a AA, Ea AA, kn AA, Ea AA kn 
AA, | Ea AA; kai AA, | Ea AA, Lu 
_ y2 i),2 (i))x(p,g)eD 
Aj = YDA) ($ ( Léi )) (p,q) 
M+1)\( N+I1 —1.2 M +l . , , 
| Jm) Pebr g Slowly converging two-dimensional 
2 2 series (similar problem arises in antenna 
N «1 
AA, = A,,,,—A, gelu : array theory [8]) 
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FSS Lattice Model 





Computation of Lattice Admittance Elements 


itta 





Procedure: 


"  |solate a single aperture from all others 


=" Solve electromagnetic aperture 
problem through classic network 
formulation [13] 


" Extract admittance for lattice model 
Original Problem Equivalence, Region A 


A B A 


E',H' E',H' 
M 


nce: single-aperture problem 


Assume equivalent magnetic 
. “ ” M VM, 
current has amplitude ("voltage") V: 


Half space self — y: — JIM, :H, (M, )dA 
admittance: À 


Input current: J, = JIM, H A4 
A 


Equivalence, Region B 





Dt TEXAS A&M 


NIV ER SI T Y 





FSS Lattice Model 


Computation of Lattice Admittance Elements 








Mutual admittance: dual-aperture problem 
"odd Equivalent magnetic current M,-V.M,, 
= Isolate two apertures from all others amplitudes V, and V;: M,-VM,, 
= Extend network formulation to two- 
aperture problem Half-space mutual — y» -||M,, -H,,(M,,)d4 
admittance: A 


= Extract admittance for lattice model 
Input currents: /,— UM. -HdA 
A 


Original Problem Equivalence, Region A Equivalence, Region B 
A B A B 
M, -M, 
i i i i 
E H EH) $ ħ 
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FSS Lattice Model 


Extension to Multiple-Element Designs 





Multiple-element FSS examples 








periodic 
i L , . 
cell ə} = Multiple resonant structures per unit 

114] cell 

Primarily achieves multiple resonant 
y bands 
l — Relatively unexplored from a synthesis 

group o - 

N apertures (or perspective 

patches) per cell By | n , ` ` 

—— Cte Stone D pe? — Can be easily cast in the lattice model 

(or patch) framework 









within group 





= Examples from literature: 
— Multiple apertures per unit cell [14] 
— Dual rectangular strips [15] 


— Multiperiodic cross dipoles (several 
layered levels of periodicity) [16] 





Reflectance 


fl 
Kl 
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FSS Lattice Model 


Extension to Multiple Element Designs 








Generic two-element example 


= Assume uniform plane wave incidence Y“ =Y + A 


— Bysymmetry, lattice collapses to two- 

rns 
port network m2k 
Two voltages corresponding to two 
apertures per unit cell 


Cr) Cr) = Two network ports per unit cell 








Y^ jer 
k 
Three different infinite summations of 


admittances 


v] [y> ye] JL, 
V,| (y> Y| |L 
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FSS Lattice Model 


Extension to Multiple Element Designs 








Generic two-element example 
| Y^! z self admittance of aperture 1 plus all 
i mutual admittances between aperture 1 
| and surrounding apertures y“ = kast Ja 


Y+? = self admittance of aperture 2 plus all y 
| 1 —V,4 sr i 
Se mutual admittances between aperture 2 m 


Cr) Cr) and surrounding apertures vi Y 


Y2? - sum of mutual admittances between k 


aperture 1 and all surrounding apertures 
of type 2 only 


lf viii 
iy^ Y>| LL 
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FSS Lattice Model 


Extension to Multilayer Designs 





A ^l. 

/ MW F F 
ac | E jJ 
Zi Ü UI 


Original Problem Eguivalence, Region A 





iittances between FSS layers 


" Two-layer design: aperture problem 
extended to three regions 
E.H E',H' — Admittances in regions A and B 
MA identical to single-laver problem 
— Two port network formulation 


= Region B between FSS screens 


— Apply image theory between two 
ground planes 


Equivalence, Region B Equivalence, Region C — Infinite number of images 
B C — Sum of images is slowly converging 


— Apply one-dimensional Shanks 
M transformation to accelerate series and 
^? 4 7 -M, estimate infinite image summation 
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FSS Lattice Model 


Extension to Multilayer Designs 





Infinite Images 
Equivalence, Region B 


" Two-layer design: aperture problem 


B 
extended to three regions 
— Admittances in regions A and B 
-M, identical to single-layer problem 
f " — Two port network formulation 
= Region B between FSS screens 
— Apply image theory between two 
4 ground planes 
— Infinite number of images 
2d d d 


— Sum of images is slowly converging 


| 

| — Apply one-dimensional Shanks 

^ sx transformation to accelerate series and 
| ! M MEN 

! estimate infinite image summation 
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FSS Lattice Model 


Extension to Multilayer Designs 





Infinite Images 
Equivalence, Region B 


" Two-layer design: aperture problem 


B 
extended to three regions 
— Admittances in regions A and B 
-M, identical to single-layer problem 
e — Two port network formulation 
= Region B between FSS screens 
— Apply image theory between two 
"e ground planes 
` ` ` ` — Infinite number of images 
A — Sum of images is slowly converging 
M | 


| | 
| | — Apply one-dimensional Shanks 
fif | 4 "tt transformation to accelerate series and 
| | . e "NP à . 
| | estimate infinite image summation 
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FSS Lattice Model 


Extension to Multilayer Designs 





Two-layer equivalent network 


" Two-layer design: aperture problem 


hs i i e 
Zu "äh "Al CH extended to three regions 
E, AY nz H AR. 


mz 


— Admittances in regions A and B 
identical to single-layer problem 





M,-V.M,, 
amplitudes V, and V;: M, -VM,, 


Equivalent magnetic current — Two port network formulation 


= Region B between FSS screens 
Half space self Y^ = -|| M., H. (M, , dA — Apply image theory between two 
A 


admittance: ground planes 
Mutual admittance Y = IM. H. (M,, )dA — Infinite number of images 
between lavers: Al — Sum of images is slowly converging 
“Self” images: mn — Apply one-dimensional Shanks 

a transformation to accelerate series and 
“Mutual” images: re » Y. (z) estimate infinite image summation 


Input current: I,- JIM, : Hi dA 
A 
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FSS Lattice Model 


Extension to Multilayer Designs 





Coupling between layers 











Region C 
"Self" image term 


Mutual coupling 
with image effects 


Region A 


Self admittance . 
Mutual coupling 
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FSS Lattice Model 


Extension to Multilayer Designs 





Collapsed lattice under uniform incidence 


Region A Region B Region C 


Y^? — self admittance of aperture 1 plus 
mutual admittance between aperture 1 
and surrounding apertures 





lou £2Y Y, 


| j=l 


Y^! z self admittance of aperture 1; 
mutual admittance between aperture 
1 and surrounding apertures; “self” 


Y^? 2 sum of all coupling admittances 
between layers (including images) 


and "mutual" image terms 
` y> = WIE, FY | 
Y^ -2Y, «2Y/ «29 Y, «Y, | ġa 
j=l 
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FSS Lattice Model 


Examples: Rectangular Aperture FSS 












| L —-7mm L,-10 mm 

















W —- 0.75 mm W=0.5mm h=4mm 
D, =D, =12 mm D, =D, =11.8mm 
8 — l g 
= -30) | — Lattice Model P. —-VI - P= 
N — . sedi d 0 | i i i 
A m Simulation ? A 2 - Lattice Model J 
— Simulation 
50 50 
T - P; T = PF T Fa? 
“Uy 15 20 25 P Sag 15 20 25 P 
Freq (GHz) [18] in Freq (GHz) in 
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FSS Lattice Model 


Examples: Rectangular Aperture FSS 





Dual layer rectangular aperture FSS 





























D. 
e—a) 
RegionA | Region B , Region C 
L=4mm | | 
>a W — 0.2 mm 
L— ——346€ W ||D, 
| D.-5 mm 
D, — 2mm 
0 | | 
p 1 HA -10 
Two identical FSS screens in jl i 
20) 
$ a g — Lattice Model 
Separation between screens: Bux y=" m a Sii 
. -4() | 
- Duroid 5880 (er = 2.2) 
T— fa 90 
- Thickness: 62 mil E. e 
710 20 30 40 50 


Freq (GHz) 
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Polynomial Representation of Admittances 
Optimization Problem 


Fabricated Examples 


FSS SYNTHESIS THROUGH LATTICE 
MODEL 
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FSS Synthesis Through Lattice Model 


Polynomial Representation of Admittances 








D "= Observations 


— Computation of lattice admittance 
L values is the most expensive part of 
— —39W | D, the model 





E eh — Lattice admittances tend to be 
a v M smoothly varying functions of the FSS 
cn NU 95 ven dimensions over a reasonable range 





" Possibilities 
d us I — Admittance variations could be 
ja Par] es. represented with good accuracy 
through multi-dimensional polynomials 
— Polynomial-based circuit models for 
specific FSS geometries 


— Efficient optimization using polynomial 
models 
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FSS Synthesis Through Lattice Model 





Polynomial Representation of Admittances 
Example of Polynomial Model 
Vary L, D,, and W simultaneously: 


Y^ (L,D,W) =cL’+c,LD+c,LW +¢,L+c.D° +c,.DW +c,D+cW° +cW +c, 


io mim Led e Ld Coefficient of determination (R2): 0.9997 
E 10.17 
D =? 
prm c, | |-187.35 Adjusted R2: 0.9996 
10 mm « D, « 20 mm C, 0.63 
0.5 mm « W «1.5 mm c; | | -133 Root-mean-sguare error - 1.354e-5 
c | | JI 
C —0.02 
Polynomial regression U 
based on 125 admittances Ce —24.29 
calculated directly C, 2.83 
Cn —0.01 
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FSS Synthesis Through Lattice Model 





Optimization Problem 


Conceptual optimization procedures 





Typical FSS optimization: 


Define multi- Nonlinear constrained optimization 

dimensional -Genetic algorithm, PSO, etc. 

optimization -Full-wave numerical computation at each 
space iteration 





'Direct' optimization using lattice model: 


Define multi- Nonlinear constrained optimization 

dimensional -Any nonlinear algorithm viable 

optimization -Full computation of lattice admittances at 
space each iteration 





"Interpolated" optimization using lattice model: 


Define multi- Compute lattice | — 
| Nonlinear constrained optimization 
dimensional admittances at small Polynomial 
ae au -Polynomial computation of admittances 
optimization number of points in regression | 
MA -Very simple implementation 
space optimization space 
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FSS Synthesis Through Lattice Model 


Optimization Problem 








Formulation of optim 


Nonlinear constrained minimization 
— Easily implemented with Matlab's 
optimization toolbox [19] 
— For simplicity, only two variables are 
varied: L and D = D, = D, 
— Transmissivity is maximized at a single 
target frequency 
Constraints 
— Linear constraints maintain realistic 
dimensions 
— Vectors x, and x, define the bounds of 
the optimization space 
— Nonlinear constraint f(x) ensures 
transmissivity over 0.95 


min -T (x) 
s.t f(x)<0 
Ax <b 
xX, SxXxXx, 


Optimization vector: 


x-[L D] 


Linear constraints: 


] -0.9 
47) | b - 
0 -l 





D 


E~ 
LL  eW||D, 


Objective function: 


P RhG) 
P, 2P,Y'(x) 


in 


7(x)= 


Nonlinear constraints: 


| f(x)20.95-T(x) 
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FSS Synthesis Through Lattice Model 


Optimization Problem 








Formulation of optimization problem: dual aperture 








= Nonlinear constrained minimax min max [-, (x).-7; (x)] 
— Easily implemented with Matlab's E f (x) <0 
optimization toolbox 
| MU Ax € b 
— Five-dimensional optimization space 
RAY xX, <x<x 
— Transmissivities at two target frequencies 
are simultaneously maximized 
" Constraints Optimization vector: Nonlinear constraints: 
e e e e e e T 
- Linear constraints maintain realistic oe WA L, D, D, n] 0.95 — T; (x) 
dimensions f (x)—| 0.95—T, (x) 
— Vectors x, and x, define the bounds of the T, (x) —0.05 


optimization space Linear constraints: 


— Nonlinear constraint f(x) ensures 


transmissivity over 0.95 at dual passbands 00 0 -09 1I -2W 
and under 0.05 between passbands A211 0 -09 0 O| b=! 0 
01 -09 0 0 0 
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FSS Synthesis Through Lattice Model 





Optimization Problem 














X-band design examples 
0 
-10 
-20 
m 
© -30 
S — Lattice Model 
ei — Simulation ` 
-50 
9 8 9 10 11 12 13 e 8 9 10 11 i2 13 
Freq (GHz) Freq (GHz) 
Target frequency: 10 GHz Target frequencies: 9 GHz and 11.3 GHz 
D | 
l L-13.26 mm L =11.55mm L, =15.56 mm 
Cer, W =1 mm W =1 mm h=16.57 mm 


D, =24.35 mm D,=16.57 mm 





D, = D, 2 26.5 mm 
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FSS Synthesis Through Lattice Model 


Fabricated Examples 





band design examples 














0 0 
-10 10 
-20; 20 
m m 
E = -30 
A0 -40 
-50 [20] -50 
9 8 9 10 11 12 13 w 8 9 10 1] 12 13 
Freq (GHz) Freq (GHz) 
" Fabrication " Measurements 
— Rogers 3850 LCP laminate substrate (2 mil — Free space fixture with X-band horns 
thickness to approximate no dielectric) — Two-port VNA measurement 
— Apertures formed by chemical wet etching — Gated-Reflect-Line (GRL) calibration [20] 





Thin FSS sandwiched between 2 mm 
Rohacell structural foam (er ~ 1.1) 


Non-trivial deviations from ideal FSS 
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Conclusion 





Summary and Future Work 


" Summary 
— New canonical element-level circuit model for FSS 
— Separates FSS analysis into two parts: 
e Electromagnetic problem of isolated unit cells 
e Circuit problem to account for periodicity 
— Straightforward extension to multi-element and multilayer designs 
— Polynomial representation of lattice admittances 
e Accurate representation of FSS over wide range of dimensional parameters 
e Simplified optimization procedure through interpolation of admittances 
" Ongoing and future work 
— Reconfigurable FSS designs: 
e Incorporation of reconfiguration mechanisms into lattice admittances 
e Optimization / controlled reconfiguration through polynomial admittance models? 
— Application to additional FSS geometries 
— Design of multi-band, multi-order filters 
— Extension of lattice model to finite FSS 


UNIVERSITY 





Dt TEXAS A&M 


M 


References 





[1] L. Fernandez-Martinez, E. Garcia, D. de Castro-Galan, and D. Segovia-Vargas, "Equivalent Circuit Model Proposal Applied to 
Frequency Selective Surfaces". 

[2] X. Yao, M. Bai, and J. Miao, "Equivalent Circuit Method for Analyzing Frequency Selective Surface With Ring Patch in 
Oblique Angles of Incidence" 

[3] R. Dubrovka, J. Vazquez, C. Parini, and D. Moore, "Equivalent circuit method for analysis and synthesis of frequency 
selective surfaces" 

[4] A. L. Gutierrez, M. Lanza, |. Barriuso, L. Valle, M. Domingo, J. R. Perez, and J. Basterrechea, "Multilayer FSS Optimizer based 
on PSO and CG-FFT" 

[5] M. Sabielny, "Design of frequency-selective radomes using parallel Particle Swarm Optimisation" 

[6] D. M. Pozar, "A relation between the active input impedance and the active element pattern of a phased array" 

[7] J. Cserti, “Application of the Lattice Green’s Function for Calculating the Resistance of Infinite Network of Resistors” 

[8] S. Hollos, “A Lattice Green’s Function Introduction” 

[9] R. McOwen, Partial Differential Equations: Methods and Applications (2" Ed.) 

[10] Y. T. Lo and S. W. Lee, Antenna Handbook, pp. 12-28 — 12-34 

[11] S. Singh and R. Singh, "On the use of Shank's transform to accelerate the summation of slowly converging series" 

[12] N. C. Albertsen, G. Jacobsen, and S. B. Sorensen, "Nonlinear Transformations for Accelerating the Convergence of M- 
Dimensional Series" 

[13] R. F. Harrington and J. R. Mautz, "A generalized network formulation for aperture problems" 

[14] J. A. Reed, "Frequency selective surfaces with multiple periodic elements" 

[15] . Delihacioglu, "Frequency selective surfaces with multiple-strip group elements" 

[16] L. Mingyun, H. Minjie, and W. Zhe, "Design of multi-band frequency selective surfaces using multi-periodicity combined 
elements" 

[17] Multilayer 

[18] Ansoft, HFSSO v15.0, Pittsburgh, PA 15219 

[19] Matlab R2012a, Math Works, Natick, MA 

[20] GRL 


Dt TEXAS A&M = nm 


UNIVERSI T Y 





